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ABSTRACT 


This  report  gives  the  user  of  the  HP41-CV  handheld 
programmable  calculator  or  the  IBM  3033  computer,  a  blade 
element  aethod  for  calculating  the  total  power  required  in 
forward,  straight  ar.d  level  high  speed  flight  for  an 
isolated  rotor.  The  computer  programs  consist  of  a  main 
program  which  calculates  the  necessary  dynamic  parameters  of 
the  main  rotor  and  several  subroutines  woich  calculate  power 
required  as  well  as  aaxiaua  forward  velocity,  stall  onset 
velocity,  and  velocity  for  best  endurance. 
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INTEOD  OCT ION 


A.  GEN3RAL 

The  basis  for  helicopter  rotor  analysis  was  developed  in 
the  early  1920's  whan  Glauert  extended  propeller  theory  to 
the  special  case  of  rotating  wings.  Since  that  tine,  the 
development  of  the  digital  computer  has  permitted  many 
improvements  on  Glauert's  analysis. 

In  order  to  develop  a  method  for  predicting  the  total 
rotor  power  required  in  forward  flight,  it  is  necessary  tc 
develop  a  method  that  accurately  predicts  rotor  dynamics. 
The  prediction  of  refer  dynamics  in  forward  flight  is  a 
complex  one.  Typically,  a  helicopter  rotor  blade  encounters 
a  flew  environment  which  changes  rapidly  as  it  moves  around 
in  azimuth. 

In  forward  flight,  rotor  blade  sections  are  subjected  to 
azimuthal  variations  in  net  only  angle  of  attach  but  also 
dach  number.  As  a  result,  comprehensive  performance 
analysis  of  a  helicopter  is  much  more  involved  than  that  of 
a  conventional  aircraft. 

Recant  helicopter  design  trends  have  been  in  the 
direction  of  increasing  the  maximum  forward  velocity 
possible  as  well  as  higher  blaie  tip  speeds.  Since  the 
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speed  of  the  helicopter  is  added  to  the  speed  of  rotation  cf 
the  advancing  blade,  the  highest  relative  velocities  occur 
at  the  tip  of  the  advancing  blade.  tfhen  the  section  Mach 
number  of  the  blade  tip  exceeds  the  critical  Mach  number  cf 
the  airfoil,  compressibility  effects  result.  These  effects 
include  a  large  increase  in  profile  drag  and  change  ir. 
pitching  moments,  therefore  creating  additional  power 
re  cuirem  ants. 

In  all  current  helicopters  there  is  a  tendency  for  the 
retreating  blade  to  stall.  Just  as  the  stall  of  an  airplane 
wing  limits  the  low  speed  perfomance  cf  an  airplane,  the 
stall  cf  a  rotor  limits  the  high  speed  potential  of  a 
helicopter.  The  relative  velocity  of  tae  retreating  blade 
decreases  as  forward  speed  increases.  However,  the 
retreating  blades  lust  produce  the  same  amount  of  lift  as 
the  advancing  blade.  Therefore,  as  the  relative  velocity  of 
j  the  the  retreating  blade  decreases  with  forward  speed,  the 

blade  angle  of  attack  must  be  increased  to  equalize  lift 
,  throughout  the  the  rotor  disc  area.  As  this  angle  increases 

j  the  blade  will  eventually  stall  at  some  forward  speed.  For 

these  reasons,  stall  and  compressibility  thus  combine  to 
limit  the  performance  fo  a  helicopter  at  higher  speeds. 
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For  efficient  design  of  a  helicopter,  the  helicopter 
designer  should  have  the  analytical  tools  necessary  tc 
predict  the  performance  of  a  helicopter.  This  report  gives 
the  user  of  the  HP4  1-CV  and  the  IBM  3033  computer  a  means  of 
calculating  the  total  power  required  in  forward,  straight 
and  level  flight. 

3.  CBJ2CTIV3 

The  objective  was  tc  provide  an  interactive  computer 
program  which  car.  be  used  during  the  initial  design  stage  to 
estimate  the  total  power  required  by  the  main  rotor  of  a 
helicopter,  as  well  as  maximum  forward  velocity  possible, 
stall  onset  velocity,  and  best  endurance  velocity. 

Additionally,  the  desired  accuracy  chosen  at  the 
beginning  of  this  effort  was  tc  obtain  angle  of  attach 


within 

a  plus  or  minus  one 

half  degree. 

power  estima 

t  e  s 

within 

ten 

per  cent  cf  the 

actual  power  r 

equired,  and 

to 

obtain 

the 

needed  accuracy 

in  as  short  a 

running  time 

as 

possible. 
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II 


PROBLEM  D  BFI NITION 


A.  DESCRIPTION  OF  PROBLEM 

For  a  helicopter  in  level  flight,  the  maximum  forward 
speed  is  limited  by  the  power  available  as  well  as  stall  and 
compressibility  affects.  It  is  therefore  advantageous  to 
develop  a  simpl e-to- use  computer  program  that  estimates 
blade  stall  and  compressibility  effects  on  helicopter 
performance. 

3.  HETHOD  CF  SOLUTION 

In  forward  flight,  the  aerodynamic  environment  of  the 
rotor  blade  varies  as  the  rotor  blade  rotates  with  respect 
to  the  direction  of  flight.  The  method  chosen  to  obtain  the 
objective  accuracy  ignores  any  variable  not  immediately 
impacting  on  the  accuracy  of  the  determination  of  angle  of 
attack  and  power  requirements. 

The  method  utilized  herein  uses  a  combination  of 
moB!en*ioi  and  blad3  element  theory  to  perform  rotor 
performance  calculations.  This  theory  is  initially  used  to 
determine  the  induced,  profile,  ar.d  parasite  power  required. 

The  solution  method  then  utilizes  a  blade  element  analysis 
to  predict  the  cyclic,  collective,  and  inflow  angles 
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associated  with  the  rotor.  These  angles  are  than  used  in 
ths  calculation  of  tha  rotor  blade’s  angle  of  attack  at  the 
azimuthal  positions  of  90  and  27  0  degrees.  Compressibility 
and  stall  power  are  then  estimated  as  a  function  of  angle  of 
attack  and  forward  velocity. 
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III.  DESCRIPTION  OF  T{J E  PROBLEM  SOLUTION 
A.  GENERAL 

The  calculaticn/pradiction  of  helicopter  perforaar.ce  is 
primarily  a  matter  of  determining  the  power  required  and 
power  available  over  the  desired  flight,  regime.  The  power 
information  may  then  be  used  to  predict  the  operational 
capabilities  of  the  aircraft. 

For  the  ideal  helicopter  (no  losses;  ,  in  forward, 
straight  and  level  flight,  power  required  can  be  subdivided 
info  five  parts:  induced  power,  PI»  required  to  produce 

refer  thrust;  profile  power,  P0  •  -squired  to  overcome  the 
sic  in  friction  and  pressure  drag  of  the  main  rotor;  parasite 
power,  Pp,  required  to  overcome  the  parasite  drag  of  the 
helicopter;  comcressibility  power,  required  to  overcome 
the  increase  in  profile  drag  when  the  tip  Sach  number 
exceeds  the  drag  divergence  Each  number  of  the  airfoil;  and 
stall  power,  Fs,  required  to  account  for  the  increase  ir. 
rotor  torque  as  a  result  of  retreating  blade  stall.  The 
HPai-cv  and  IBB  3033  computer  programs  contained  herein 
provides  a  simple,  quick  means  of  evaluating  the  power 
required  by  an  isolated  rctor. 
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3.  ASSUMPTIONS 


The  major  task  in  helicopter  performance  analysis  is  the 
determination  cf  rotor  forces,  angles,  and  power.  The 
method  chosen  ignores  any  variable  that  does  not  directly 
impact,  on  the  desired  accuracy.  Therefore,  in  order  to 
simplify  the  computational  process  used,  the  following 
assumptions  are  made: 

1.  steady  flow  through  the  rotor  system. 

2.  Small  angle  approximations  are  a  valid 
representation  of  the  real  world  phenomena. 

3.  All  blades  considered  are  rectangular  (non- tape  red) 
with  cnly  uniform  twist  possible. 

4.  Hinge  off-set  is  zero  (i.e. ,  me  thrust  vector 
passes  through  the  C.G.). 

5.  The  stall  angle  or.  most  helicopter  blades  can  be 
approximated  by  the  angle  that  occurs  at  CLMAX  for  a  2-D 
airfoil. 

6.  The  helicopter  is  trimmed.  This  implies  that  the 
sum  of  all  the  moments  about  the  center  of  gravity  (C.G.)  is 
zero,  all  forces  are  in  balance,  and  that  no  lateral 
flapping  is  present. 
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7.  At  stall  onset,  :hr  value  3 f  sac-ion  drag 
coefficient  jumps  approximately  0.08. 

8.  The  rotor  rotates  counterclockwise. 

C.  NOMENCLATURE 

"Standard”  nomenclature  is  used.  Appendix  A  contains  an 
alphabetical  list  defining  all  the  symbols  and  mnemonics 
used  in  the  development  of  this  report.  Appendix  3  contains 
an  alphabetical  list  of  all  HP* 1-cv  displays  used  in  the 
HPhl-CV  program. 

D.  INITIALIZATION 

It  is  assumed  that  the  initial  design  of  the  helicopter 
has  been  completed  and  that  the  helicopters  wsigat  as  well 
as  the  chord,  radius,  ti?  velocity,  twist,  zero-lift  drag 
coefficient,  and  number  of  blades  of  tae  rotor  are  known. 
Finally,  it  is  assumed  that  an  initial  estimation  of 
equivalent  fiat  plate  area  is  available,  and  than  the 
forward  velocity  of  the  helicopter  is  known. 

E.  ROTO 8  DYNAMICS 

There  are  four  rotor  parameters  which  will  help  expedite 
later  calculations.  The  first  of  these  is  the  rotor  advance 
ratio,  a  ;  u  is  the  ratio  of  the  helicopter's  forward 
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(2) 


u  *  V, 


T 

w  a  ar  a  , 


*  (  n  R) 


The  second  dime r.sionless  ratio  that  needs  to  be 
calculated  is  the  inflow  ratio,  \  .  The  inflow  ratio  is  the 
ratio  cf  the  net  velocity  up  through  the  rotor  system  to  the 
tip  speed.  ?cr  the  near  hover  case,  u  <  3.1,  the  Tan  a3  3 
and  the  inflow  ratio  can  be  approximated  by 


(3) 


/ctTT 


In  forward  flight,  the  calculation  of  X  requires  the 
determination  cf  the  angle  of  attack  of  the  rotor  disc, 
letting  the  angle  of  attack  between  -he  iisc  plans  and  the 


-*! 


incoming  free-stream  velocity  be  a3  ,  and  assuming  angles  to 
be  small,  the  inflow  velocity  ran  be  caloulated  using 

(3a)  X  =  -  cT  ♦  u  Tan  a3 

2/X2+y2 

w  h  er  e 

a3  =  -Tan  (Dp/L  >  =  -Tar.  (0p/V) 

Dp  =  <Pp*  550)/  7f 

The  last  dimensionless  parameter  naan  needs  to  be 
calculated  is  solidity.  Solidity  is  tns  fraction  of  the 
rotor  disc  area  that  is  composed  of  blades.  For  a  blade  of 
constant  chord  (i.a.  ,  non-tapered)  solidity  can  be  expressed 

33 

(4)  J  *  (b  *  C)/(PI  *  S) 

Finally,  the  tip  loss  factor,  3TL,  *us”  135  considered. 
The  tip  less  factor  is  used  to  account  for  the  loss  of  lift 
that  a  rotor  blade  experiences  due  to  flow  from  the  rotor’s 
lower  surface  to  its  upper  surface.  The  tip  loss  factor  of 
a  rotor  can  be  approximated  by 

(5)  Btl*  1.0  -  A  *  CT  /  b 
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P.  VELOCITY  CALCULATIONS 


There  are  four  velocities  that  are  of  interest.  The 
first  of  these  is  the  downwash  velocity  ,  w.  Assaminq 
steady  flow  through  the  rotor,  the  downwash  velocity  can  be 
approximated  by 


(6)  v 


a 


2P  An  V, 


(ft/s) 


NOTE:  This  equation  is  not  valid  for  small 
values  of  forward  velocity. 


The  second  velocity  that  needs  to  be  calculated  is  the 
stall  onset  velocity,  Vg  .  The  stall  onset  velocity  is  the 
velocity  at  which  the  retreating  blade  tip  first  exceeds  the 
static  stall  angle.  The  forward  speed  at  which  stall  onset 
is  first  noted  can  be  approximated  by  toe  velocity  for  best 
range,  VBR.  This  is  due  to  the  marked  increase  in  profile 
power  required  at  speeds  higher  than  velocity  for  best 
range.  A  typical  set  of  power  curves  for  a  helicopter  are 
shown  in  Figure  3.1. 

The  forward  velocity  for  minimum  ?/’/  (i.  a.,  best  range 
velocity)  is  easily  found  graphically  on  the  power  required 
curve  as  the  point  where  a  straight  line  tn rough  the  origin 
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POWER  REQUIRED  VS  VELOCITY 


Pigure  3.7.  Typical  Helicopter  Power  Curves 


is  tangent  to  the  curve  (see  Figure  3.2).  Since  th?  pove 
curves  are  r.ot  initally  known,  an  analytical  expression 
needs  to  be  developed  which  will  estiaate  the  velocity  f c 
best  range.  The  following  is  the  derivation  of  cts  such 
expression. 
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As  shown  below,  if  the  parasite  power  required  is  set 
equal  to  the  profile  power  required  and  the  equality  is  then 
solved  for  an  equation  in  which  forward  velocity  is  the 
variable,  the  result  is  a  cubic  equation  in  which  the 
largest  rcot  defines  the  point  where  me  profile  power  and 
parasite  power  are  equal. 

D  x  V 

■  P  -0 

1/2Pfy=  1/8  a  Cio  p  Ad(1  «•  4.25  y*  ) 

(7)  y-  (4.25  /4 f)  Cdc  Ad7t  -  (0/4  i)  cio  kDy  *  o 

Next,  the  largest  root  of  equation  7  needs  to  be 
determined.  This  can  be  simply  done  using  the  cubic  roct 
equation.  Letting 

p  =  -  (4.25  o/4f)  Cio  Ad7t 
r  *  -  (  a  /4  f )  cdo  advt3 

Equation  7  can  now  be  written  as 
(7a)  y  -  ?«*  -  r  *  0 
Substituting 
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(  X  -  p/3  ) 


for  Vp  in  equation  7a  yields  ar.  equation  of  the  fora 

(7b)  X3  ♦  aX  ♦  b  =  0 

where , 

a  =  -1/3  ?2 
b  =  1/27  (2p*  ♦  27r) 

The  largest  root  of  equation  7b  is  then,  qiven  by 

(8)  V  *  A  ♦  3 

w  h  ar  e , 

A 

3 

Equation  8  can  now  be  used  as  the  initial  apcrcxiaat ion 
stall  onset  velocity. 

The  third  velocity  which  is  of  interest  is  ‘he 
helicopter's  maxiaua  forward  velocity,  /MAX.  As  shown 


=  [-+  *  ] 
*  [=5-  -  V**  '*-*♦  ] 


1/3 


1/3 
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Assuming  that  -he  variation  of  profile  power  with  forward 
velocity  i s  negligible,  the  velccity  for  best  endurance 
(also  the  best  rate  of  climb  velocity)  ran  be  found.  If 
equation  is  differentiated  with  respect  to  forward 
velccity,  VF,  and  is  set  equal  to  zero,  it  can  be  seen  that 

(11  a)  ?x  *  3  Pp 

or 

(11b)  '42  =  3  0  f  Vp2 

2  p  ADVBE  *■ 


Solving  equation  for  the  best  endurance  velocity,  7 BF, 
r e  su its  in 


(11c) 


(ft/s) 


5.  INITIAL  POWER  CALCULATIONS 

As  forward  velocity  increases,  the  induced  power 
decreases,  the  profile  power  increases  slightly,  and  the 
parasite  power  increases  until  it  becomes  the  dominant  loss 
at  high  speeds  [Ref.  2].  For  forward,  straight  and  level 
flight,  the  induced  power  can  be  calculated  by 
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l 


(12) 


If  tip  losses  are  taken  into  affect  (aquation  5), 
induced  power  now  becomes 


(12a)  ?j  =  (  1/3^)  *  ?x 
(TL) 

Since  the  induced  power  required  in  ground  effect  i 
than  that  required  cut  of  ground  effect,  equation  12 
be  written  as 


(12b)  Pz  *  (1/%,)  *  (G2)  *  ?! 

(TL+GE) 
where. 


G2  =  (-0.  1276  *  (h/D)  *  +  0.708  *  (h/D)  3 
-1.4569  *  (h/D)  2  ♦  1.3432  *  (h/D) 

♦  0.5  147) 


The  profile  power  required  is  given  by 

(13)  ?0=aCdoPAD^  (1  ♦  4.25  U2> 

4400 


Finally,  parasite  power  given  by 
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50 

(h?) 

the 


L  r  SS 

should 


(hp) 


(14) 


?p  •  Jf  ¥ 

1100  (hp) 

H.  ANGLE  CF  ATTACK 

The  determination  of  angls  of  attar*  at  the  azimuthal 
positions  of  90  and  270  degrees  are  important  in  the 
determination  of  compressibility  at i  stall  effects.  It  is 
because  of  this  effect  that  the  dynamics  of  the  blade  motion 
are  important  in  analyzing  helicopter  performance. 

The  angle  of  attack  of  the  rotor  is  a  function  of 
radius,  r,  and  azimuthal  position,  ip  .  Figure  3.3 
illustrates  the  sign  convention  used  in  determining  blade 
angle  of  attack.  The  angle  of  attack  cf  a  rotor  can  thus  be 
estimated  by. 


(15) 

a  (r ,  \p  ) 

0  -  VpBCos a  ♦  *  - 

Va3 

ft  R  ♦ 

vFSin  * 

(r: 

■*  h  er  e , 

(15a) 

0 

=  V 

0T+  0^CO3  'F 

♦  Osin  'F  ♦  K  S 

2  O 

(15b) 

6 

=  3o" 

al  Cos  i|;  - 

bl  Simp 
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ANGLE  OF  ATTACK 


To  determine  angle  of  attach,  it  is  taan  necessary  to 

determine  the  longitudinal  collective  and  cyclrc  ingles 

and  0  .  This  can  be  acconolished  bv  expressing  the 

2 

coefficient  of  thrust  as 
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(16)  C-  *  i£  ATI  ♦  (0  ♦  KnS)  T2  «-0_T3  ♦  (8  -  K_b1)T4 

i  2  0  O  X  2  D 

wher ? , 

T1  «  .5  (32  ♦  .  5  u 2>  T3  *  .  259*< 3 2  ♦  u2  ) 

TL  TL  TL 

T2  *  (.  5  B3  ♦  .  5  u2  3  )  T4  *  .5  y  (32  ♦  .25  u  2) 

TL  TL  TL 

Additionally,  the  longitudinal  flapping  coefficient  al, 
needs  to  be  determined.  The  longitudinal  flapping 
coefficient  can  be  written  as 


(17)  al  a  A  All  ♦  (  0  ♦  A12  +  9^13  ♦  (8  -  K  b1)11« 


where , 


D1  =  (B3  -  .5  U3> 
TL 


All  a  .25  (.5  ^2l-  y3  /3) 

B~  2  D  1 
TL 


A 1 3  -  2  y  B  2 

M  TL 


DT 


Al  2  =  8  V  3, 


TL 


3  D  1 


A14  =  Bt£  *  1.5u2 

- jH - 


Assuming  that  there  is  nc  lateral  flapping  and  that  the 
:ect  is  zero  (i.e.,  KB  *  0.3),  agnations  16  and  17  become 
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(18) 


2  C<v  =  ATI  ♦©  P2  ♦  9  T3  +  0  T4 


(19)  al 


X  All  ♦  0  A  1 2  ♦  3™A  1  3 

o  * 


+  3  A14 

2 


NOTE : 


Since  the  analysis  is  only  longing  at 
the  azimuthal  ocsitions  of  90  ar.d  270 
decrees,  it  car.  be  seen  iro.ii  equation 
15b  that  the  ccr.tr  i  ban  ion  of  tne 
longitudinal  flapping  coefficient,  al, 
wilx  always  be  zero  at  these  cositions. 


Equations  18  and  19  new  represent  a  set  of  si multanecus 
equations  in  which  the  only  unknowns  are  the  collective  and 
cyclic  angles  and  can  thus  be  determined.  Knowing  the 
values  of  the  cyclic  and  collective  angles,  the  blade  tip 
angle  of  attack  at  the  90  and  270  degree  positions  can  be 
estimated  by 


(20) 


(21) 


„  = 

0  + 

0m  +  0  +  X 

90 

0 

T  2  1+U 

x~-n  = 

0  - 

0  +  0  +  A 

2  /0 

0 

^  T  — 

NOTE: 

The 

DO  SI 
UD  , 

aperie  of  attack, 
tively  rf  the  li 
see  Figure  3.2. 

;se 
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I.  STALL  POWSH 

In  the  helicopter  stall  normally  starts  at  the  tip  of 
the  retreating  blade,  since  the  highest  angles  of  attack  are 
usually  at  the  blade  tip.  As  the  forward  speed  increases, 
the  stalled  area  of  the  rotor  blade  spreads  inboard. 

At  the  higher  values  of  v,  the  effects  of  stall  or.  newer 
required  are  great  and  therefore  need  to  be  estimated. 
Assuming  a  jump  of  0.08  in  the  value  of  Tin  at  stall  onse*, 
tnat  the  rotor  area  within  which  the  blade  stall  exists  is  a 
segment  of  minimum  dimensionless  radius  Xg  and  that  the 
stall  area  is  symmetric  abcut  <P  *  270  degrees,  Castles  and 
New  found  that  the  effects  of  tip  stall  on  power  required  at 
the  higher  values  of  u  are  large  and  can  be  approximated  for 
high  speed  flight  by 

<22>  cps  3  -snr  <’  -1“  >*  <’  -  y  /rrT 

where  Xs  is  the  non dimensional  radius  outboard  of  whicn  the 
retreating  blade  is  stalled  (Ref.  3].  The  dimensionless 
radius,  Xg  ,  can  be  estimated  by  equating  the  section  angle 
of  attack,  at  P  =  270  degrees,  to  Ad  AX  *Ref.  «].  Setting 
equation  15  at  p  -  270  degrees  equal  to  AMAX  results  in  the 
quadratic  listed  below. 
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(23) 


(X  -u) 


Qo(x  -y)  ♦  3t  f-  Q^y-e^x  -u)  ♦  >.*x  ji 


Applying  -the  quadratic  foraula  to  equation  23  yields  the 
roots 

(23a)  Xs 

(23  b)  X 

0 

(23c) 

(23b) 

(23  s) 

Equation  22  is  satisfactory  for  nest  cases.  It  is 
possible  however,  for  the  blade  section  angle  of  attack  to 
be  higher  inboard  than  at  the  tip  creating  a  situation  which 
is  ussuaily  refered  to  as  inboard  stalling  [Bef.  5]. 

For  the  special  case  of  inboard  stalling,  see  Figure 
3.4,  the  incremental  stall  power  coefficient  defined  by 
equation  22  is  too  large  and  needs  to  be  corrected. 
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STALLING 


t 


INBOARD 


Figure  3.4.  inboard  Stall  Patern. 


.\ssuming  a 

he  stalled  region 

is 

diamond  3 .taped  as 

shown  in 

Figure  3.5 

and  that  the  stall 

ed 

area  is  symmetric 

about  'h 

3*/2,  it 

can  be  seen  that 

as 

X  aoDroaches  X  , 
o'  S 

*  V>  15 

correction 

to  the  incremental 

st 

all  power  coeffici 

CPS 

must  vanish. 

34 


STALL  REGIONS 


Si miliarly ,  as  the  average  of  X  and  aooroaches  a 

o  S  ■  * 

"ha  value  or  Cpg  sees  to  the  value  defined  by  equation 
Therefore,  in  order  to  correct  for  the  possibiii 
inboard  staliina,  the  correction  *a c*o--  k  .  :  ? 


(24) 


1.0 


k  = 
S 


1.  ) 


(25)  k  =  B  /2G_  ♦ 

S  T  S 

■^1  "-T£  " 


1.0 


The  corrected  equation  for  stall  newer  required  thus  becomes 


(26) 


C* 

PS 


*  c 


PS 


J.  COMPR3SSI3ILITY  P0W35 


The  individual  affects  cf  stall  and  compressibility  on 
rotor  profile  power  are  substantial  at  nigh  advance  ratios. 
When  both  effects  are  present  the  losses  lua  to  each  source 


i  ire  difficult  to 

di.  3 

anguish.  Therefore, 

as  a  helicopter’s 

\  forward  velocity 

and 

tip  speed  increase. 

the  need  for  a 

simple  estimate  of  how  the  com pressibilr tv  of  the  air 
influences  the  rotor  performance  is  necoessary. 

In  forward  flight,  the  Mach  number  cf  the  advancing 
blade  is  given  by 
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(27) 


a  =  sTIp  (x  ♦  u  sir.  *  ) 


where , 


(27a) 


(27b) 


m  *  ns  =  •/ 

TIP  T 


x  =  r/R 


nondimensional  radius 


Since  t he  highest  Mach  number  occurs  at  one  tip  of 
advancing  blade  at  y  -  90  decrees,  equation  27  ca: 


M  ( 1  ♦  u  ) 
TIP 


Reference  3.5  shewed  that  the 
divergence  car.  be  estimated  by 


ML  =  MCRO  -  ml  *  a  *  a 
CRIT 


Mach  number  for  dr' 


Gesssow  and  Cria, 


cheir  investigation 


compressibility  found  tha*  the  compressibility  effect  on 
rotor  performance  was  a  rapid  increase  in  the  profile  pewe: 
when  tha  tip  Mach  number  exceeded  the  critical  Mac.o  number 
for  drag  iivergen.es  [Ref.  6].  The  increase  in  profile  pews: 
coefficient  due  to  Mach  effects  car.  be  estimated  oy 
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CDMCLiJS  ICSS  AND  3  ZCOflil  EN  D  ATION  S 


Th~  objective  of  this  project  was  to  provide  an  easy  -o 
use,  inter-active  computer  program  which  could  be  used 
during  the  initial  design  phase  zo  estimate  the  coral  power 
required  by  the  main  rotor  of  a  helicopter.  The  c:a?u*e: 
programs  contained  herein  provide  results  wnich  are  well 
within  the  objective  accuracies  (see  Figures  C.l,  Z.Z,  1.31 
and  provide  acceptable  results  as  a  first  cut  estimate  of 
compressibility  and  stall  power  requirements. 

In  the  development  of  the  computational  model,  many 
simplifying  assumptions  were  made  to  ease  the  amount  of 
computation  required.  The  assumption  which  most  impacts  on 
the  accuracy  of  the  program  is  that  of  steady  flow  *hrough 
the  rotor. 

The  flow  environment  encountered  by  the  rotor  changes 
rapidly  Jue  to  the  rate  of  change  cf  clads  angle  cf  attach.. 
Additionally,  rotor  operation  at  high  advance  ratios  also 
produces  considerable  radial  flow  along  the  blade  span  . 
The  steady  flow  assumption  ignores  the  time-variant  aspect 
of  rotor  aerodynamics.  The  dynamic  nature  of  ta a  rotor, 
especially  when  operating  at  or  near  me  stall  regime, 
requires  the  application  of  unsteady  aerodynamics,  and  a 
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close  examination  of  how  the  pitching  moments  get 


:h«  recreating  blade  stalling  affects 
controiability  [Bef.  7].  It  is  therefore  recomended  " 
additional  investigations  consider  how  the  unsteady 
aerodynamics,  and  pitching  moments  generate!,  influence 
performance  of  the  helicopter. 
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APPENDIX  A 
MCMEHCLA  TUP.3 


Term 

Mnemonic 

De  f irit ion 

a 

a 

Slope  of  airfoil 

a  action  lift  curve. 

a 

AD 

5  of  or  disc  araa. 

a3 

Aipha3 

Disk  plane  angle  of 

attack. 

aMAX 

AMAX 

The  steady  flow  stall 

angla  of  the  airfoil 

(given  by  a  *  CL  MAX  ] 

a 

90 

A90 

Angle  of  attack  cf 

the  advancing  blade 

at  <|>  «  90  degrees. 

a 

270 

A270 

Angle  of  attack  cf 

the  retreating  blade 

at  'li  ~  27 0  deg  rees. 

al 

al 

Longitudinal  flapping 

All 

AM 

Term  in  definition  of 

THETA2. 

M2 

A  12 

Term  in  definition  of 

Units 


rad 


rad 


rad 


rad 


rad 


rad 

dimensionless 


dimensionless 
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I 


HETA2 


A  1  3 

A  14 

b 

BET  AO 

B5 

8TL 

C 

CDO 

C  L  M  AX 

CPC 


'PS 


ws 


CPS 


cs 


Tara  in  definitor.  of 
THETA2 . 

Tarn  in  definition  of 
THETA2. 

Number  of  blades  on 
the  main  rotor . 

Sain  rotor  coning 
angle. 

Ter  a  in  definition  of 
XS  * 

Tip  loss  factor. 

Sean  chord  of  main 
rotor  blade. 

Sain  rotor  coefficient 
of  drag  at  zero  lift. 
Saxiaua  cceficier.t  of 
lift  (2  -  D)  . 

Correction  to  powar 
coefficient  due  to 
compressibility  effects. 
Correction  to  power 
coefficient  due  to  stall. 
Tara  in  definition  of  Xg. 
42 


1 


dimensionless 


drae  r.sionless 


dimensionless 


rad 


dimensionless 


dimensionless 


dimensionless 


dime  nsior.iess 


diaensionlsss 


dimensionless 


dimensionless 


I 


CT 


oefficiert  of  thrust. 


aimer. 


C 


T 


Dp  D 

D  DP 

AS-  DSD 


5  DEL  3 

3 


f  F?  A 


i 


i 


GA!1MA 


G.  E.  G 
h  H 


hp  H? 

K  KB 


of  the  main  rater. 

Bator  iisc  diameter. 
Parasite  drag  of  the 
helicopter. 

Term  in  the  definition 
of  compressibility 
?  o  w  er . 

Sate  of  change  of  blade 
pitch  with  respect  to 
blade  flapping. 
Equivalent  flat  plate 
area  of  the  helicopter 
ir.  forward  flight. 

Term  in  definition  of 
3  an  d  C  . 

Ground  effect  ratio. 
Height  of  main  rotor 
above  the  ground. 
Horsepower, 
affect. 


lb 


aimer 


Limer 


iimen 


aimer 


hp 


dimer 


X 


LAMB 


M  MACH 

M  __  MCRIT 

CRIT 


M  MT 

T 


’^CRTO 


MCHTO 


MU 


ml  Ml 


Q  0 8 SGA 

?  PM 

C 


■5 

‘  I 


PI 


Ratio  cf  the  net 
velocity  up  through 
the  rotor  system  to 
the  tip  speed. 

Mach  number  of  rotor 
blade. 

Critical  Mach  number 
of  advancing  blade 
at  <p  »  90  degrees. 
Tip  Mach  number  of 
rotor  balda. 

Critical  Mach  number 
for  Cl  =0.0  . 

Main  rotor  advance 


dimensionless 


dimensionless 


dimensionless 


dimensionless 


dimensionless 


dimensionless 


rat  io. 

Constant  in  definition 


of  critical  Mach  number 

(ml  =  0.  113)  . 

Rotational  velocity. 

rad/s 

Pover  required  due  to 

hp 

compressibility  effects. 

Induced  power. 

**  P 
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•o 

o 

?0 

Profile  power. 

hp 

0 

'  P 

P? 

Parasite  Power. 

hp 

?s 

PS 

Power  required  due  fc 

hp 

stall  effects. 

R 

R 

Main  rotor  radius. 

P 

RHO 

Air  density. 

slug/ft 

a 

SO 

Main  rotor  solidity. 

listens! 

0 

THETA 

aatio  of  ambient  temp¬ 

dimetsi 

erature  to  standard  sea 

level  temperature. 

0 

THETA3 

Main  rotor  collective 

rad 

0 

pitch. 

Main  rotor  lateral 

rad 

pitch. 

0 

THETA2 

■Main  rotor  longitudinal 

rad 

2 

cyclic  pitch. 

T1 

T 1 

Term  in  definition  of 

a ime  ns : 

THETAO . 

T2 

T  2 

Term  in  definition  of 

dimers: 

THETAO. 

T3 

T  3 

Term  in  definition  of 

diasns: 

THETAO. 
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T4 

Term  in  definition  of 

dimers: 
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t 


r  H  E  TA  0 


VI 

71 

Induced  velocity. 

ft/s 

7f 

7? 

Aircraft  forward  speed. 

ft/5 

V 

7  XT 

Aircraft  forward  speed. 

ht 

VM  AX 

VMAX 

Maximum  forward  velocity 

possible. 

Kt 

7s 

V  3 

Snail  onset  velocity 

(velocity  at  which  A270 

equal  A  MAX) 

hr 

VT 

V  T 

Rotor  tip  speed. 

ft/3 

a 

SV 

SVEL 

Soeed  of  sound. 

ft/3 

w 

DW 

aotor  dcwnwash  velocity. 

ft/3 

W 

W 

Aircraft  gross  weight. 

lb 

xs 

X3 

Radius  outboard  of  which 

the  Bain  rotor  is  stalled. 

dimensionless 

X 

0 

xo 

Radius  inboard  of  which 

rotor  blade  stall  aav 

be  present  due  to  inflow 

ratio  and  blade  twist. 

dime  nsionless 

a6 


APPENDIX  3 

SI  i lz£l  ?scga.\a  documentation 

A.  GENERAL 

The  HP  41-CV  program  uses  45  program  storage  registers. 
Prior  to  program  initialization ,  the  handheld  calculator 
should  be  sized  to  46. 

The  computer  program  consists  of  a  main  program  which 
calculates  the  necessary  dynamic  parameters  of  one  main 
rotor  and  several  subroutines  which  calculate  power  required 
as  well  as  the  maximum  forward  velocity,  stall  onset 
velocity,  and  best  endurance  velocity. 

The  documentation  for  the  HP  4  1-CV  program  is  deviled 
into  the  following  sections: 

1.  PURPOSE 

This  section  describes  the  intended  purpose  of 
the  program  or  subroutine. 

2.  ASSUMPTIONS 

This  section  lists  an y  assumptions  made  which  are 
applicable  to  the  program  cr  subroutine. 

2.  EQUATIONS 

This  section  lists  the  equations  utilized  within 
the  mam  orogram  cr  subroutine.  The  primary  references  for 


the  aquations  used  are  Aerod vnanics  and  V/STO L  Flight. 
Reference  3.1,  deli,  cop  ter  Th  e  or  y .  Reference  3.2,  and 
Aircraft  Performanc a.  Reference  3.3. 

4.  FLOWCHART 

Both  the  handheld  computer  and  the  IBM  3033 
no  mail-/  execute  instructions  in  a  program  in  a  sequential 
manner  unless  it  is  instructed  to  do  otherwise.  This 
section  will  graphically  represent  the  step  oy  step  method 
used  to  solve  the  problem  as  well  as  the  flow  of  control 
between  the  various  parts  of  the  program.  In  a  flowchart, 
different  types  cf  operations  are  indicated  by  different 
shaped  boxes  as  illustrated  below: 


Oval 


Recta  ngl e 


For  start  or  stop. 


For  a  calculation  c: 
process  other  than  • 
decision. 


5.  PROGRAMS  AN D  SUBROUTINES 


This  section  lists 

all  the 

pro  grams 

subrcuti 

nes  which  must  be  present 

in  the  HP 

41-cv  ?rio 

program 

a  xecution. 

6. 

PROGRAM  LISTING 

This  section  contains  t h 3  3?  41-cv  listing  of  the 
program  or  subroutine. 

B.  QUICK  REFERENCE  TA3LES 

The  tables  in  this  section  will  0?  useful  tc  the  user  of 
the  H?  41-cv  handheld  calculator  as  a  source  of  quick 
reference.  Table  1  is  an  alphabetical  listing  of  all 
calculator  displays  with  an  explanation  of  their  respective 
meaning.  Table  2  lists  all  storage  registers  used  and 
describes  how  they  are  utilized. 
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< 


1 


•  H?41  -CV  Disolavs 

TABLE  t 


or 

SPLAY 

definition 

7NTT 

a 

=  ? 

Slops  of  airfoil  lift  car  vs. 

rad 

b 

=  ? 

Number  of  main  rotor  blades. 

:i/a 

c 

-  7 

Chord  length  of  main  rotor. 

DA 

S  ? 

Density  altitude. 

c  • 

Cl 

c  =  ? 

Sain  rotor  coefficient  of  drag 

N /A 

at  zero  lift. 

F? 

A  =  ? 

Equivalent  horizontal  flat  plate 

^  — 

area  of  the  helicopter. 

h 

=  7 

Height  of  the  main  rotor  aoova 

6  ± 

the  ground. 

P  0  =? 

Critical  Mach  number  for  the 

N/A 

coefficient  of  lift  equal  to 

zero. 

3 

-  7 

Main  rotor  radius. 

TW 

1ST  =  ? 

Geometric  twist  of  the  rotor. 

rad 

7? 

(KT)  =  ? 

Forward  velocity. 

xt 

7T 

-  ? 

P.ctor  tip  speed. 

ft/ 

4 

-  7 

Aircraft  gross  weig.it. 

Los 

CD 

T  *  — 

Main  rotor  ccllecriy3  tit  oh. 

ra  3 
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blade  at  1  =  270  degrees. 


2.  H?ui -cv  Register  Utilization 

TABLE  II 


STORAGE 


QISTER 

STORED  j 

00 

Ground  e£ 

01 

Maxiaum  f< 

02 

Stall  ons< 

05 

Sot  or  rad. 

06 

Number  of 

07 

Zero- lift 

09 

Rotor  h<?i' 

10 

Air cr  aft 

1  1 

Air  den  si 

12 

Lift  curv 

13 

Rotor  tip 

14 

Coerf icie 

15 

Tip-loss 

16 

.lair,  rcro 

17 

Rotor  hei 

19 

Coaipr  assi 

i 


o  (h/d)  . 


i 
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19  Solidify  (SD). 

20  Induced  velocity  (VI). 

21  Main  rotor  profile  power  (P0). 

22  Advance  ratio  (MO)  . 

23  Sr  a  11  power  (PS)  . 

24  Advance  ratio  squared  (MU2  )  . 

25  Forward  velocity  (  VF)  . 

2b  Equivalent  horizontal  flat  plate  area  (??.\) 

2~>  Maximum  2-D  lift  coefficient  (CLdAX)  . 

28  Main  rotor  parasite  power  ( P P)  . 

29  Main  rotor  geometric  twist  (TWIST). 

ui  Angle  of  attack  at  27  3  degrees  (A2?3). 

42  Angle  of  attack  at  90  degrees  (A93)  . 

43  Sonic  velocity  (SV2L). 

44  Critical  Mach  number  for  coefficient  o?  lif 

equal  to  zero  (MCRO). 

45  Density  altitude. 

OTHERS  Scratch  pad  calculations. 
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C.  PROGRAM  documentation 

This  section  contains  the  necessary  documentation  :or 
the  H?  41-CV  computer  program.  The  main  program  as  well  as 
all  of  the  subroutines  used  in  the  solution  of  the  problem 
are  outlined  in  “his  section. 
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1.  3 Program 

(a)  PURPOSE 

This  program  calculates  the  dynamic  parameters  of 
the  main  rotor  which  are  necessary  for  calculating  the 
total  main  rotor  power  required  in  forward,  high  speed 
straight  and  level  flight.  It  additionally  controls  the 
execution  sequence  of  the  various  subroutines  which  are 
used  to  calculate  main  rotcr  power  required,  in  terms 
of  horsepower,  as  well  as  maximum  forward  velocity,  stall 
onset  velocity  and  velocity  for  best  endurance  in  icnots. 

(b)  ASSUMPTIONS 

(1)  All  angles  are  small. 

(2)  Steady  flow  through  the  rotor. 

(3)  All  rotor  blades  are  rectangular  (non-tapered) 
with  only  uniform  twist  being  possible. 

(4)  Only  the  first  hamonir  of  flapping  is  necessary 
for  calculating  power  required. 

(U)  The  effective  dimensionless  radius  can  be 
apprcxima teay  by  the  tip-loss  factor. 

(6)  The  thrust  vector  passes  through  the  C.G.. 

(7)  The  static  stall  angle  for  blades  on  most 

helicopters  can  be  approximated  by  the  angle  at 

which  CLNAX  occurs  for  the  2-D  airfonl. 
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(8)  The  <$  3  affect,  or  the  result  of  cockir. 
flapping  axis  of  the  blade  so  that  its  ? 
varies  as  the  blade  flaps  is  zero  (i.e. , 
lateral  flapping  is  present). 


)  EQUATIONS 


7  *  1.68394 


'4  *  VVT 

VF  (ft/S)  * 

T1  *  .5(Bt£  ♦  .5  u*> 

T  2  =  .  3  34B  3  ♦  .5  u*B 

TL  TL 

T3  =  -25Bt£<  EJy*  ♦  y*> 

T  4  =  .  5  U  (8  *  ♦  .25  y2) 

TL 

D1  =  <3t]?  -  .5  u*> 

All  =  4  (.5  Bt£  -  u 3  /3)  /D 1 

A  12  *  8uBtl/3D1 

A 1 3  *  2u^*  /  01 

A 14  *  (B  *  ♦  1.5  u«>  /  D1 
TL 

( 2CT/aa)  =  ATI  ♦  9^12  ♦  9^,13  +  9^4 

0  *  XA11  ♦  9  A12  ♦  0  A 1 3  +  9  A14 

o  T  2 

Anax  *  CLNAX/CLA 


a27J  =9  -  9 

0  2 

F  =  AMAX  -0 


cs  = 


ur  +  x 

-  u9t-  r 


x.  ■  -3„» 


s 

X  * 


<V 


■3s-  <V 


♦  9  ♦  (  X  / 1  ♦  U  ) 
-0 


-  49TCs)-3  /  29t 

-  49tCs)-s  /  29t 
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PROGRAMS  AND  SUBROUTINES  USED 


f.  PROGRAM  LISTING 


01*LBL  ‘MBS’ 
02  •?='?* 

03  PROMPT 
04  ST0  05 
05  -b=7' 

06  PROMPT 
07  STO  06 
08  -C=7* 

W  PROMPT 
10  STO  04 
!1  -Cd0=7* 

12  PROMPT 

13  STO  07 

14  -W^- 

15  PROMPT 

16  STO  10 
! 7  -VT=?- 

18  PROMPT 

19  STO  13 
20  ■vf<KT>='> 

21  PROMPT 

22  1.68894 

23  * 

24  STO  25 

25  *FPfl=,>" 

26  PROMPT 

27  STO  26 

28  •«=?* 

29  PROMPT 

30  STO  09 

31  ’a*’* 

32  PROMPT 

33  STO  12 
34  -CLMflX=7’ 

35  PROMPT 

36  STO  27 
37  -TWISTS'?* 

38  PROMPT 

39  STO  29 
*«  *MCR0=7* 

41  PROMPT 

42  STO  44 


43  -Jft=7* 

44  PROMPT 

45  STO  45 

464LBL  -PGH" 

47  RCL  25 

48  RCL  13 

49  / 

50  STO  22 

51  X92 

52  STO  24 

53  xec  -sd- 

54  XEQ  -BEH* 

55  FS7  03 

56  GTO  05 

57  XE8  -PIT- 

58  XEQ  -POl- 

59*LBL  05 
60  XEQ  -PP1* 

61  RDV 

62  FS7  03 

63  GTO  06 
64  XEQ  -VBE- 

65  XE9  -VM0X- 
66  QDV 

67*LBL  06 

68  XEQ  -CT- 

69  RCL  15 

70  YT2 

71  STO  30 

72  ENTERt 

73  RCL  24 

74  2 

75  /' 

76  * 

77  STO  31 

78  2 

79  / 

30  STO  33 
81  RCL  15 


32  3 

83  YtX 

84  3 

85  / 

86  EHTEPt 

87  RCL  2* 

38  RCL  15 

89  * 

90  2 

91  / 

92  ♦ 

93  STO  34 

94  RCL  38 

95  RCL  24 

96  + 

97  RCL  30 

98  * 

QQ  4 
100  / 

101  STO  35 
182  RCL  38 

103  RCL  24 

104  4 

105  / 

106  + 

107  RCL  22 
188  * 

109  2 

110  / 

111  STO  36 

112  RCL  22 

113  RCL  30 

114  * 

115  2 

116  / 

117  RCL  22 

118  3 

119  YTX 

120  8 
121  / 

122  - 
123  4 


61 


16b  RCL  12 

287  - 

124  * 

167  / 

288  / 

125  RCL  38 

168  RCL  19 

289  STO  32 

126  EHTERt 

169  / 

218  FS7  03 

127  RCL  24 

178  STO  42 

211  GTO  07 

128  2 

171  XEQ  -LRffi- 

212  R-9 

129  / 

172  RCL  03 

213  -CYCLIC: 

138  - 

173  RCL  33 

214  «RCL  X 

131  STO  32 

174  * 

215  RVIEU 

132  RCL  30 

175  CHS 

216  STOP 

133  * 

176  RCL  42 

134  t 

177  ♦ 

217*LBL  87 

135  STO  37 

178  RCL  29 

218  RCL  31 

136  RCL  22 

179  RCL  35 

219  EHTERt 

137  RCL  15 

188  * 

228  RCL  36 

138  * 

181  - 

221  RCL  32 

139  8 

182  STO  31 

222  * 

140  * 

183  RCL  83 

223  - 

141  3 

184  CHS 

224  RCL  34 

142  / 

185  RCL  37 

225  / 

143  RCL  32 

186  * 

226  STO  33 

144  / 

187  RCL  29 

227  FS’  83 

145  STO  38 

188  RCL  39 

228  GTO  88 

146  RCL  22 

189  * 

229  R-D 

147  2 

190  - 

230  *COLL=’ 

148  • 

191  STO  32 

231  RRCL  X 

149  RCL  30 

192  RCL  34 

232  RVIEM 

150  * 

193  * 

233  STOP 

151  RCL  32 

194  EHTERt 

234  QBV 

152  ' 

195  RCL  31 

153  STO  39 

196  RCL  38 

2354101  08 

154  RCL  30 

197  * 

236  RCL  33 

155  EHTERt 

198  - 

237  RCL  32 

156  RCL  24 

199  RCL  34 

238  - 

157  1.5 

200  EHTERt 

239  RCL  29 

158  * 

281  RCL  40 

240  ♦ 

159  ♦ 

282  * 

241  EHTERt 

160  RCL  32 

283  RCL  36 

242  RCL  83 

161  / 

284  EHTERt 

243  1 

162  STO  40 

285  RCL  38 

244  RCL  22 

163  RCL  14 

286  * 

245  + 

164  2 

165  * 
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246 

/ 

247 

•f 

248 

STO  41 

249 

PS?  03 

250 

CTO  09 

251 

R-D 

252 

*<270=' 

253 

PRCL  X 

254 

PVIEH 

255 

STOP 

256 

PIV 

257  XE9  -1%’ 

258  PBV 

259  RCL  27 
288  RCL  12 
261  / 

262  RCL  33 

263  - 

264  RCL  32 

265  + 

266  STO  30 

267  RCL  22 

268  CHS 

269  RCL  29 

270  * 

271  RCL  38 

272  - 

273  STO  31 

274  RCL  30 

275  RCL  22 

276  * 

277  RCL  03 

278  + 

279  STO  35 

280  4 

281  CHS 

282  * 

283  RCL  29 

284  * 

285  STO  35 

286  RCL  31 


287 

298 

289 

290 

291 

292 

293 

294 

295 
2% 

297 

298 

299 

300 

301 
382 

303 

304 

305 

306 

307 
388 


Xt2 

327UBL  83 

n  I  l. 

STO  36 

328  XE8  ‘CPS' 

EHTERt 

329  OBV 

RCL  35 

330  CF  81 

CHS 

331  CF  82 

x>r 

332  CF  83 

SF  01 

333  XEC  *CPC 

FS?  81 

334  PDV 

GTO  83 

335  PDV 

RCL  36 

336  XE9  -PT! 

RCL  35 

337  XEQ  ‘CHG 

+ 

SORT 

338*LBL  89 

STO  38 

339  XEO  *VS 

RCL  31 

340  EHB 

CHS 

+ 

RCL  29 

/ 

2 


STO  30 


309  ’XS=" 

310  PRCL  X 

311  RVIEU 

312  STOP 

313  RCL  31 

314  CHS 

315  RCL  38 

316  - 

317  RCL  29 


318  / 

319  2 

320  / 

321  STO  88 

322  *XO=* 

323  PRCL  X 

324  RVIEU 

325  STOP 

326  PBV 
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2.  Soli it* v 


8.  PURPOSE 

Thts  subroutine  calculates  the  ratio  of  total  blade  area 
to  the  total  rotor  disc  area. 

0.  ASSUMPTIONS 
None 

C.  EQUATIONS 

b  *  c  *  R  ft  *  c 

°  '  pi  *  R* ~  ’ 


d.  FLOWCHART 

^start 


calculate 
sol Idlty 

r'  ■ 

^return 
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PftOGHAMS  AND  SU8RCUUNES 

•  W8S  * 

1.  PrtOGAAM  L  ISMNG 


i 


0HLBI  -SP- 
02  BCL  06 
03  RCL  M 

04  * 

05  RCl  05 
06  ' 

07  PT 

09  3T0  !9 
16  EHB 


USED 
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a.  PURPOSE 

The  purpose  of  this  subroutine  is  to  compute  the  induced 
velocity  of  a  rotor  system. 

b.  ASSUMPTIONS 

Steady  flow  through  the  rotor  system. 

c.  EQUATIONS 

*  - 

2<>A0  v, 

d.  FLOWCHART 


e.  PROGRAMS  ANO  SUBROUTINES  USED 

"  was* 
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PROGRAM  LISTING 


B1*L8L  -W 
02  RCL  25 
03  X=0? 

04  GIG  05 
05  RCL  10 
06  2 

0?  RCL  1! 
00  * 

09  RCL  25 
10  * 

11  RI 

12  * 

S3  RCL  05 

14  Xt2 

15  * 

16  t 

17  STO  64 

18  GTO  06 

miBL  05 
20  SF  03 

2!»LBL  06 
22  EHl 
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4. 


Coeificii^  of  Thrust 


a. 


PURPOSE 

This  subroutine  calculates  tne 
arbitrary  rotor. 


coefficient  of  thrust  for 


o.  ASSUMPTIONS 

Steady  flow  through  me  rotor  system, 
c.  EQUATIONS 


d.  flowchart 


C  start 

■  i  “ 


_ _ i - - 

calculate  j 

coeff  of 
thrust  1 
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e.  PROGRAMS  AND  5UBPOU TINES  USED 


•  was  • 

f.  PROGRAM  LISTING 


81*L8L  *CT* 
02  9CL  85 
@3  Xt2 
04  PI 
05  * 

06  PCL  11 
07  * 

08  RCL  13 
0?  Xt2 
10  * 

11  1/X 

12  PCI  10 

13  * 

14  STO  14 

15  END 
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5.  Induced  Power 

(a)  PURPOSE 

This  subroutine  calculates  the  power  required  by  the  rc 
to  produce  thrust  at. hover  and  forward  flight.  Additionally, 
this  subroutine  corrects  for  tip  losses  (losses  in  lift  at  the 
tips  due  tip  votices)  as  well  as  ground  effect. 

(b)  ASSUMPTIONS 

Steady  flow  through  the  rotor  systea. 

(C)  EQUATIONS 


-  1.4569(h/D>  *+»34  3  2(h/0|  +0.5147 
PIT  =  ( 1/B  )*(OE)  *PI 
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e.  PROGRAMS  ANO  5U0ROUI INES  USED 


*  was  ‘ 

f.  PROGRAM  LISTING 


8 1 ♦LBL  -PIT- 
02  RCL  14 
03  2 

@4  * 

05  SWT 
00  RCL  00 
07  / 

08  CHS 
09  1 

11  STO  15 

12  RCL  09 

13  RCL  05 

14  2 

15  * 

10  ' 

17  STO  17 

18  1.5 

19  XOY 

20  x>r 

21  GTO  95 

22  RCL  17 

23  1.3432 

24  * 

25  RCL  17 

26  Xt2 

27  -1.4509 

28  * 

29  + 

30  RCL  17 

31  3 
3?  YtX 


33  .  788 

34  * 

35  * 

36  RCL  17 

37  4 

38  Y4X. 

39  -. 1270 

40  * 

41  + 

42  .5147 

43  ♦ 

44  STO  90 

45  CTO  96 

404LBL  95 
4?  1 

48  STO  99 

494LBL  00 

50  RCL  25 

51  Xt2 

52  STO  32 

53  RCL  10 

54  2 

55 

56  RCL  11 

57  / 

58  B! 

59  / 

60  RCL  95 

61  Xt2 

62  ' 


63  SWT 

64  STO  20 

65  Xt2 

66  1/X 

67  RCL  32 

68  * 

69  2 
79  / 

71  CHS 

72  STO  34 

73  Xt2 

74  l 

75  * 

76  SWT 

77  RCL  34 

78  + 

79  SWT 
88  RCL  10 
81  * 

82  RCL  20 

83  * 

84  550 

85  / 

86  RCL  15 

87  * 

88  RCL  09 

89  * 

90  STO  16 

91  -PI=- 

92  9RCL  X 

93  RVIEH 

94  STOP 

95  END 
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6.  Profile  22112 

a.  PURPOSE 

This  subroutine  calculates  the  prof  tie  poser  required  for 
forward,  strafgnt  and  level  flight  tn  terms  of  horsepower. 

O.  ASSUMPTIONS 

Steady  flow  through  the  rotor  sustem. 


C.  EQUATIONS 


P„  *  o  Cdo  p  A  V*  (1  ♦  4.25  u:) 

0  0  T 


d.  flowchart 


4400 


/" 


start 


’  . 

calct. 

prol 

PO|l 

ilate 

'Me 

»er 

, 

r 

store  PO 

Tn  R21 

(np) 


return 


e.  PROGRAMS  ANO  SUBROUTINE  USED 
•  WBS  • 
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PROGRAM  LISTING 


01»LBL  *P01 
02  RCL  19 
03  RCL  07 
04  * 

05  RCL  11 
0b  * 

07  RCL  05 
00  yt2 

09  * 

!0  oi 

11  * 

12  RCL  13 

13  3 

14  Y+X 

15  * 

16  8 


18  STO  21 

19  RCL  24 

20  4.25 

21  * 

22  1 

23  ♦ 

24  RCL  21 

25  * 

2b  550 

27  / 

28  STO  21 

bQ  -P0r‘ 

30  0RCL  X 

31  0VIEU 


7.  Parasite  Power 


a.  PUPOSE 

This  subroutine  calculates  tne  parasite  power  required  In 
forward,  straight  and  level  flight. 

b.  ASSUMPTIONS 

Steady  flow  through  the  rotor  system. 

c.  EQUATIONS. 


PP 

d.  FLOWCHART 


»fV* 


1100 


Chip) 


start 


_ 

calculate 

parasite 

power 

store  PP 

In  R28 

, 

t 

return 


«.  PROGRAMS  ANO  SUBROUTINES  USED 
•  W8S  " 


f5 


T 


PROGRAM  LIST  TNG 


81*LBL  'Bp!’ 

92  PCI  11 

93  PCL  26 

94  * 

95  .5 
86  * 

97  »CL  25 

98  3 
09  YtX 
10  * 

11  550 

12  / 

12  3TQ  28 

14  *PPs* 

15  9RCL  X 

16  9VIEH 
1?  STOP 
18  ENB 
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8.  Maximum  ?cr ward  V ^locl* y 


a.  Purpose 

This  subroutine  calculates  the  power-1 Imlted  maximum  speed 
of  the  specified  he if copter. 

t>.  ASSUMPTIONS 

a.  The  power-limited  maximum  velocity  may  be  estimated 
by  neglecting  the  varlatln  of  Induced  power  and  profile 
power  with  speed. 

b.  Power  required  to  hover  Is  approximately  equal  to 
power  required  fo  maximum  speed. 

c.  Steady  flow  through  the  rotor  system. 

c.  EQUATIONS 

VMAX  -  V,  tJL  ]  'n 


V 


f 


2  P 


V 


1/2 
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flowchart 


PftOUxAMS  AND  SUBROUTINES  USED 


•  was  • 

PROGRAM  l 1ST  INC 

81*LBL  *VRW 

82  PCI  26 

83  SI 

*4  / 

85  PCI  05 
8b  Xt2 
8?  / 

88  1'X 

89  4 
■8  * 

11  .3333334 

12  VtX 

13  PCL  20 

14  * 

15  $70  81 
lb  1.6884 

17  / 

18  *V«flX=- 

19  PRCL  X 

20  PVIEii 

21  STOP 

22  PCL  81 

23  PCL  25 

24  x>r 

25  GTO  12 

26  X<=*7 

27  GTO  13 

284LBL  12 
-vf  >  VflftX" 
30  PVIEW 
31  XE9  *CNG- 

32*LBL  13 
33  END 
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9.  Bes*  Hr. iura  n.ce  V  elcci  ~v 


A.  PURPOSE 

To  calculate  the  value  of  velocity  corresponding  to 
minimum  power  (t.e.,  best  endurance  velocity  and/or  best 
rate  of  climb). 

B.  ASSUMPTIONS 

1.  Steady  flow  through  the  rotor  system. 

2.  The  variation  of  profile  power  with  forward  velocity 

is  negligible. 


C.  EQUATIONS 


0.  FLOWCHART 


..  r  W  .  _ _ & —  i/i 

VK  *  C  p  *  A  *  (3  *  f  }  3 


start 


_ 

calculate 

VBE 

1 

convert 
VBE  to 
knots 

, 

return 


(ft/s) 
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e. 


PROGRAMS  AND  SUBROUTINES  USED 


•  WBS  ‘ 

f.  PROGRAM  LISTING 

*V8E‘ 

'.;fjCL  85 

,  7 

H4  ?! 

05  * 

96  3 
07  / 

08  ?CL  26 
00  / 

10  SORT 

11  RCL  10 

12  * 

13  ?CL  *1 

14  ■ 

15  RCL  05 

16  m 

i” 

18  P1 
1?  " 

20  SORT 

21  1.6880* 
22  ' 

23  *VBE*=* 

24  RRCL  X 

25  QVIEW 

26  STOP 

27  END 


01 


10. 


Stall  Onset  Velo city 


a.  PURPOSE 

This  subroutine  gives  the  user  of  the  HP  41-CV  an  Initial 
approximation  for  the  velocity  at  which  the  retreating  blade 
angle  of  attack  Is  approximately  equal  to  the  Is  approximately 
equal  to  the  static  stall  angle  of  the  rotor  blade. 

b.  ASSUMPTIONS 

1.  Steady  flow  through  the  rotor  system. 

2.  Stall  onset  velocity  Is  approximately  equal  to  the 
velocity  for  best  range  (t.e.,  minimum  P/V). 

c.  EQUATIONS 

p  =  (-4.25  eCdo  A  VT  )/(4f ) 

O  T 

r  -  (-7  Cdo  Ao  Vj)/(4f) 

a  =  -1/3  *  if 

b  =  1/27  *  (2  *  p’«-  27  *  r) 

A  -  t  -b/2  +  (b’/2  ♦  aJ/27)l'*l  1/3 

B  -  C  -b/2  -  (b */2  ♦  aJ/27)1/i]  1/3 


32 


d.  flowchart 


e.  PROGRAMS  ANO  SUBROUTINES  USED 
•  WBS  " 
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f.  PROGRAM  LISTING 


81*tBL  ‘VS* 

82  FS9  83 

83  GTO  18 

84  RCL  18 

85  RCL  87 

86  * 

87  PI 

88  * 

89  RCL  85 

18  XY2 
11  * 

12  4 

13  / 

14  RCL  26 

15  / 

16  STO  83 

17  RCL  13 

18  * 

19  4.25 
28  * 

21  CHS 

22  STO  81 

23  RCL  13 

24  3 

25  Yt* 

26  RCL  03 

27  * 

28  CHS 

29  STO  02 
38  27 

31  * 

32  RCl  01 

33  3 

34  Vtjr 


38  27 

39  ' 

40  STO  31 

41  RCL  01 

42  X.t2 

43  -3 

44  10! 

45  * 

46  3T0  38 


47  RCL  31 

48  X*2 

49  4 

58  / 

51  RCL  38 

52  3 

53  YtX 

54  27 

55  / 

56  4 

57  SORT 

58  STO  33 

59  RCL  31 

60  CHS 

61  2 
62  / 

63  ♦ 

64  . 333334 

65  YTX 

66  STO  32 

67  RCL  31 

68  CHS 

69  2 

70  / 

71  RCL  33 

72  - 

73  .333334 

74  YtX 

75  STO  34 

76  RCL  32 

77  + 

78  STO  82 

79  ,25 

80  * 

8!  RCL  02 
82  + 

S3  STO  25 
84  SF  03 
85  GTO  'WH* 

864LBL  10 

87  RCL  27 

88  RCL  12 
99  ' 

90  STO  31 

91  RCL  41 

92  - 


93  XT2 
9<  SORT 

95  .801 

96  XOY 

97  X<=Y9 

98  GTO  11 

99  RCL  31 
180  RCL  41 
10!  X>Y’ 

102  GTO  12 

103  RCL  41 

104  RCL  31 

105  X>Y9 

106  GTO  13 

1074L8L  12 

108  RCL  25 

109  5 

110  - 

111  STO  25 
112  GTO  -PGN- 

1134LBL  13 

114  RCL  25 

115  5 

116  * 

117  STO  25 
118  GTO  -«G0- 

1194LBL  11 
128  CF  93 

121  RCL  41 

122  P-P 

123  *4'STBLL=- 

124  ORCL  X 

125  PVIEM 

126  STOP 

127  RCL  25 

128  1.68894 

129  / 

130  -VS=- 

131  ORCL  X 

132  PVIEM 

133  STOP 

134  EMU 


11.  inflow  Rat  lo 

(a)  PURPOSE 

This  subroutine  calculates  the  ratio  of  on  a  r.er  velo 
through  the  rotor  system  tc  the  tip  speed. 

(b)  ASSUMPTIONS 

Steady  flew  through  the  rotor  system. 


(c)  EQUATIONS 

X  =  -  /CT/2 

Xs"  ^  ♦  u  Tar.  a 

‘  v'  A"  +  U‘  3 

Dp  *  ( ?P*  5  50)  /V p 

a  3  =  -Tan  (Ep/»> 
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FLOWCHART 


return 


e.  PROGRAMS  AND  SUBROUTINES  USED 


•  WBS  ‘ 

f.  PROGRAM  LISTING 

♦LBL  -LtWB' 

*2  RCL  1* 

03  2 
0*  / 

05  SORT 
96  CHS 
0?  STQ  03 
00  .1 
00  RCL  22 
10 1 «:<=r 

11  GTO  00 

12  RCL  28 

13  558 

14  « 

15  PCL  25 

[i  . 

17  RCL  10 

18  ■ 

19  STO  23 


20*lSL  0! 

21  RCL  14 

22  CHS 

22  RCL  03 
2*  *72 

25  RCL  24 

26  ♦ 

27  SORT 

28  2 

29  * 

30  / 

31  ENTERt 

32  RCL  23 

33  R-B 

34  TflH 

35  RCL  22 

36  * 

37  ~ 

38  STO  08 


39  RCL  03 

40  - 

41  Xt2 

42  SORT 

43  .00005 

44 

45  GTO  02 
*6  RCL  08 

47  STO  03 

48  GTO  01 

494LBL  02 
58  RCL  08 
5!  STO  03 

524LBL  06 
53  EHB 
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.** 


12.  Angle  c f  Annack  at  9  0  Deaii£3 


a .  PURPOSE 

This  subroutine  calculates  the  angle  of  attack  at  the 
azimuthal  position  of  90  degrees. 

b.  ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

b.  Blade  oscillations  are  periodic  In  nature. 

c.  Only  first  harmonics  of  flapping  are  necessary  for 
calculating  angle  of  attack. 

d.  The  thrust  vector  passes  through  the  C.G. 

e.  Only  uniform  twist  of  the  rotor  blade  Is  possible. 


C.  EQUATIONS 


□90  =  e  ♦  e  ♦  e  ♦  ~r~.  r 

O  2  t  (  l  ♦  y  ) 


i 


return 


PROGRAMS  ANO  SUBROUTINES  USED 


•  W8S  * 
"  a90  " 


PROGRAM  LISTING 


gi»LBL  *490* 
82  RCL  33 
03  PCL  32 

04  t 

05  PCL  29 
06  * 

07  EWTERt 
0S  RCL  03 
09  ! 

10  »CL  22 

11  ♦ 


14  2TP  42 

*15  R-fl 

16  * 4905s* 

17  OREL  X 

18  aVIEW 

19  STOP 

20  EHB 
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13.  Com 


or ?ssibiXit v  Fow*r 

a.  PURPOSE 

This  subroutine  calculates  the  power  required  due  to  com_ 
presslbllfty  on  the  main  rotor  sustem  In  forward, straight 
and  level  flight  In  terms  of  horsepower. 

2.  ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

b.  The  compress  1 bll Ity  losses  can  be  expressed  as  a  function 
of  the  amount  by  which  the  drag  divergence  Mach  nunber  Is 
exceeded  at  the  tip  of  the  advancing  blade. 

C.  EQUATIONS 


AM  -  M.  -  M  .  0.06 

d  t  crlt 

M  »  0.71  -  2.3  *  «90 

erf  x, 

C  -  *  CO. 012  *  AM  ♦  0.1  *  (AM  )  *1 

pc  d  d 

P  •  C  *p*PI*tf*V*  (HP) 
x  pc  w  t 

550 
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f.  PROGRAM  LISTING 


81*lBL  ‘CPC’ 
02  RCL  13 

83  RCL  25 

84  + 

85  RCL  13 

86 

0?  STO  31 
88  RCL  12 
08  RCL  12 
10  * 

11  .113 

12  * 

13  CHS 

14  RCL  11 

15  + 

16  STO  32 

17  CHS 

18  RCL  31 

to  * 

20  .86 

21  * 

22  STO  33 

23  8 

24  Vf'Y 

25  ¥<=VP 

26  GTO  09 

27  PCI  33 

28  3 

29  ytx 

30  .1 

31  * 

32  £HTERf 
S3  RCL  33 
31  .012 
35  « 

7  A  *■ 


37  RCL  1? 

38  * 

79  cto  31 

10  RCL  1! 
4!  * 

42  RCL  05 
13  X*2 

44  * 

45  R! 

46  * 

47  ENTER* 

48  RCL  13 

49  7 

50  W 

51  * 

52  558 

53  / 

54  STO  18 

55  GTO  10 


56*LBL  0? 

57  8 

58  STO  18 

59*LBL  10 

60  -pC=* 

61  QPCL  X 

62  QWIEW 

63  STOP 

64  EHD 
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14. 


Stall  Power 


a .  PURPOSE 

This  subroutine  estimates  the  additional  power  required  In 
forward,  straight  and  level  flight  due  to  retreating  blade 
stall.  Additionally,  this  subroutine  calculates  a  stall 
correction  factor,  k  that  corrects  for  the  special  case, 
of  Inboard  sta 1 1 Ing. 


b.  ASSUMPTIONS 

a.  Steady  flow  through  the  rotor  system. 

b.  The  section  drag  coefficient  at  stall  jumps  approxi¬ 
mately  0.08  at  stall  onset. 

c.  The  stalled  area  Is  symmetric  about  the  270  degree 
azimuthal  position. 

d.  For  all  airfoils  considered,  the  static  stall  angle  Is 
approximately  12.5  degrees. 


c.  EQUATIONS 


-B 


2  *  Q 


s  1.0 


where. 


B  «  -u  *  e  -  T 

•  i 

T  -  0.218166  -  e  ♦  e 

o  i 


•<*»-(  JL 


B  /2  *  9  ♦  X 


T  ■ 


1  -  X 


c  .  c  2T~Tp7*'  {1  -  "  >  *  a  -  n  *  (1 


2  .Vi 

.  '  ] 
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P  =  C  *  p  *  PI  *  R*  *  V* 

So*  T 


3li 


I  AO- A 1 36  376  COMPUTER  PROGRAMS  FOR  HELICOPTER  HIGH  SPEED  FLIGHT 
ANALYSIS(U)  NAVAL  POSTGRADUATE  SCHOOL  MONTEREY  CA 
W  F  CARMONA  SEP  83 


UNCLASSIFIED 


F/G  1/3 


NL 


MICROCOPY  RESOLUTION  TEST  CHART 

NATIONAL  BUREAU  OF  STANDARDS- 1963-A 


FLOWCHART 


PROGRAMS  ANO  SUBROUTINES  USED 
•  WBS  * 


f.  PROGRAM  LISTING 

01*LBL  ’CPS* 

02  PS?  01 
03  GTO  00 
04  RCL  30 
05  RCL  08 
06  + 

07  2 
08  / 

08  1 

10  XOY 

11  X<=Y'? 

12  GTO  01 

13  1 

14  STO  23 

15  GTO  02 

16*LBL  01 

17  RCL  31 

18  2 

19  / 

20  RCL  29 

21  / 

22  RCL  30 

23  + 

24  CHS 

25  EHTERt 

26  1 

27  RCL  30 

28  - 

29  f 

30  STO  23 


31M.BL  02 

32  "KS=* 

33  RRCL  X 

34  AVIEH 

35  STOP 

36  SF  02 

37  RCL  19 

38  2* 

39  / 

40  PI 

41  / 

42  EHTERt 

43  1 

44  RCL  22 

45  - 

46  Xt2 

47  • 

48  1 

49  RCL  30 

50  - 

51  * 

52  1 

53  RCL  30 

54  XT2 

55  - 

56  SORT 

57  * 

58  RCL  23 

59  * 

60  STO  23 

61  RCL  13 

62  3 

63  YtX 

64  RCL  23 

65  * 

66  PI 

67  * 

68  RCL  05 

69  XtZ 

70  * 

71  RCL  11 


72  * 

73  550 

74  / 

75  STO  23 

76  FS?  02 

77  GTO  07 

78»LBL  06 

79  0 

80  STO  23 

81M.BL  07 

82  -PS=’ 

83  mi  X 

84  RVIEW 

85  STOP 

86  END 
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15.  Total  Power 


a.  PURPOSE 

This  subroutine  calculates  the  total  power  required  In 
forward,  straight  and  level  flight,  to  Include  stall  and 
compressibility  power  In  terms  of  hosepower. 

C.  ASSUMPTIONS 

Power  losses,  such  as  transmission  and  cooling,  can  be 
Ignored. 

C.  EQUATIONS. 

P  *  P  *  P  +  P  +  P  ♦  p 

T  f  •  p  9  M 

d.  FLOWCHART 
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6.  PROGRAMS  AND  SUBROUTINES  USED 
*  WAS  * 

f.  PROGRAM  LISTING 


01MJL  -PTi- 
02  RCL  28 
03  RCL  16 
04  + 

05  RCL  21 
06  ♦ 

07  RCL  23 
08  + 

W  RCL  18 
10  ♦ 

11  *PT=* 

12  flRCL  X 

13  flVIEU 

14  STOP 

15  END 
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16.  Der.si.y/Sor.ic  Vslcci^v 


a.  PURPOSE 

To  generate  values  of  air  density  and  sonic  velocity. 

b.  ASSUMPTIONS 

(1)  Geopotentlal  altitude  (H),  and  geometric  altitudes  are 
equal  below  20,000  feet  (actual  ah  *  29  ft). 

(2)  In  the  troposphere  the  standard  temperature  lapse 
rate  Is  -3.57  V  per  1000  feet. 


C.  EOUATIONS 

e  -  T/T  »  (1.0  -  6.8753  E-06  *  H  ) 

SSL 

SVEL  *  SVE^  *  SQRT(e) 

RHO  -  0.0023769  *  (1.0  ♦  HTH  *  (-.02875  ♦  0.000275  *  HTH) 
where, 

HTH  »  DA/1000 


d.  FLOWCHART 


^return 
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e.  FROGMANS  AND  SUBROUTINES  USED 


*  WBS  * 

f.  PROGRAM  LISTING 


81*L8L  -BEN- 
82  RCL  45 
83  6.875  E-86 

84  * 

85  CHS 

86  1 

87  + 

88  SORT 

89  1116.89 
18  * 

11  STO  43 

12  RCL  45 

13  1888 

14  / 

15  STO  33 

16  .888275 

17  * 

18  -.82875 

19  ♦ 

28  RCL  33 
21  * 

22  1 

23  ♦ 

24  .8823769 

25  * 

26  STO  11 

27  END 


TOO 


t7.  Char,  as 


a.  PURPOSE 

This  subroutine  Is  used  to  expedite  the  changing  of  up 
to  five  of  the  Input  parameters  whenever  a  design  restraint 
Is  exceeded  and/or  at  the  end  of  the  main  program. 

b.  ASSUMPTIONS 

The  primary  parameters  which  will  require  changing  are: 


a.  Forward  velocity,  VF;  (Kts) 

b.  blade  twist,  TW;  (rads) 

c.  lift  curve  slope,  a;  (per  rad) 

d.  tip  velocity,  VTj  (ft/sec) 

e.  weight,  (lbs) 


C.  EQUATIONS 

VF  (ft/sec)  »  VF  (KtS)  *  1.68894 
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d.  FLOWCHART 


yes 


set  flags 


04 


f  lap 


DISPLAY: 


VF(Kts) 


prompt : 
VF(Kts) 


cnange  Kts  to 
ft/sec  c  store 
_ B2S _ 


start 


PROGRAMS  AND  SUBROUTINES  USED 


•  MBS  * 

PROGRAM  LISTING 


81*L8L  -CHG- 
82  'AGN* 

03  flSTO  31 

04  CF  84 

05*LBL  06 
06  ‘CHANGE’- 
07  PROMPT 

88  X=0? 

89  GTO  07 

10  SF  04 

11  SF  27 

12  -VF  TH  a  VT  H- 
13  PROMPT 

14*LBL  A 
15  *VF<KTS>a’* 

16  PROHPT 

17  1.68894 

18  * 

19  STO  25 

20  GTO  85 

21»LBL  B 
22  -THIST=’* 

23  PROMPT 

24  STO  29 

25  GTO  05 

26*LBL  C 

27  -a=’- 

28  PROHPT 

29  STO  60 

30  GTO  05 
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31*LBL  B 

32  *VT=’* 

33  PROHPT 

34  STO  13 

35  GTO  85 

36»LBL  E 

37  •«=?■ 

38  PROHPT 

39  STO  10 

40H.BL  05 

41  CF  27 

42  GTO  06 

43*LBL  07 
*4  FS?  04 
45  GTO  IMS  31 
46  END 


I 


D.  HP  41-CV  SAMPLE  OUTPUT 

Demanding  on  the  value  of  forward  velocity  there  are 
three  possible  types  of  output  that  the  HP  4 1-cv  computer 
prograa  is  capable  of  producing.  Examples  of  the  three 
possible  cases  are  listed  below. 


CASE  1  :  Forward  velocity,  V?,  is  less  than  stall 
onset  velocity. 


♦ 


PI»79i.751l53 
P0=224. 584453 
PP=8. 888888 

V$= 137. 762848 
VBEK=?8. 868750 
VHHX=i63.443338 

CYCLIC=0.888888 
C0LL=17. 896907 

4278=4.388790 

i98D=4. 888798 

P?=8. 888888 

PN=8. 808808 


PT=1  ,.016.255686 
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( 


C&SE  2 


:  Forward  velocity,  VF,  -s  3 

stall  onset  velocity. 


c. 


ter 


•ha 


?I=113.87 

P0=352.43 

PP=704.24 

VS=137.76 

VBEK=76.06 

VMflX=183.45 


CYCLIC *-9.39 
C&L=20.76 


*270=17.75 


49011s”  1.03 


*8=0.74 

X0=1.24 

*8=0.96 

PS=140.90 

PM253.71 


PT«l,573.16 


1 05 


CASE  3 


is  great? 


\ 

f 

i 


t 


Forward  velocity,  VF, 
the  maximum  forward  velo 


PM07.18 

P0=368.92 

f>P=844.?2 

VS=!3'.76 
VBEk=76.06 
Vm=l62.*5 
vp  >  mx 

CYCLIC=~!9.55 

COLL=22.08 

4278=19.87 

49811=- 1.23 

<S=8.73 

*0=1.53 


<5=1.00 

PS=!56.18 

PW=289.23 


PT=i . 766. 13 


1  06 


» 


ify  po 


APPENDIX  C 

A.  IBM  3033  PROGRAM  DOCUMENT  ATI  ON 

This  program  calculates  the  power  required  fly 
helicopter  in  forward,  straight  and  level  high  speed  fiigh 
In  the  sections  mat  follow  all  inputs  and  implement ati 
requirements  are  specified.  The  problem  solving  n e* he d ole 
used  is  as  described  in  Chapter  II. 
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3.  INPUT  DATA  REQUIRED 

Ir.  calculating  the  performance  of  a  helicopter  in  is 
necessary  to  define  a  set  of  force  conditions,  environmental 
conditions,  and  physical  conditions.  The  force  conditions 
required  are  the  aircraft  gross  weight  (W)  in  pounds,  the 
maximum  lift  coefficient  (CL.1AX)  ,  and  the  the  coefficient  of 
list  at  zero  angle  of  attack  (CDO)  .  The  environmental 
conditons  include,  forward  velocity  (VFK»  in  knots,  speed  of 
sound  (SVEL)  in  feet/sec,  rotor  height  above  one  oround  (H) 
in  feet,  and  air  density  (RHO)  in  lb-ft  /sec  .  For  sake  of 
user  simplicity,  the  the  spaed  of  sound  and  air  density  are 
generated  within  the  program  and  are  comparable  to  the 
values  found  in  a  standard  atmosphere  table.  Finally,  the 
Dhysical  conditions  required  are  the  rotor  radius  (S)  in 
feet,  tip  velocity  (VT)  in  ft/sec,  nuaoer  of  blades  (b)  , 
main  rotor  chord  (C)  in  feet,  flat  plate  area  (FPA)  in  ft  , 
geometric  twist  of  the  rotor  (TWIST)  in  radians,  and  airfoil 
lisft  curve  slope  (CLA)  in  per  radians. 
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IBM  PROGRAM  FLOWCHART 


1  10 


f 


0 


calculate 
paras 1  te 


calculate 


calculate 

A11.A12, 

A13.A14 


calculate 

Initial 


calculate 

Inflow 


calculate 
cyclic, 
coi lectlve 


calculate 
stal  1 

parameters 


CPS  «  0 

calculate 

XO.XS 

calculate 

CPS 


calculate 

KS 


correct  1 on 


i 

r  1 

KS 

-  1 

calculate 

a90 


2.  IBM  PROGRAM  LISTING 

This  section  contains  tha  listing  of  the  I3M 
computer  program  developed  in  this  report. 
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SAMPLE  OF  IBM  COMPOTES  OUTPUT 


This  section  contains  an  axample  run  of  the  IBM 
program#  utilizing  AH1-J  Cobra  fata,  starting  at  a 
velocity  cf  120  knots  and  terminating  at  VMAX. 
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oomput 
f  orwa 


i 


UGH  SPEED  FQRmARD  FLIGHT  ANALYSIS 


INPUT  DATA 


RAOIUS  •  • . . . 

MAIN  ROTOR  CHGRC... . 

NUMBER  OP  MAIN  POT OR  BLADES . 

A I  CRAP  T  GRCSS  w  EIGHT  . . 

ROTOR  TIP  VELOCITY . 

HORIZONTAL  PLAT  PLATE  AREA . 

NUMBER  OF  ENGINES  IN  HELICCPTER.. 

DENSITY  ALTITUDE . . 

MAIN  ROTOR  HEIGHT  ABOVE  GROUND... 

AIR  OENSITY  (RHCI . 

SUNIC  VELOCITY.  . . 

BLADE  GEOMETRIC  TWIST . 

TYPE  AIRFOIL . . . 

MAXIMUM  2-D  LIFT  COEFFICIENT..... 
.2tQ-.ST  ATI.C  STALL  ANGLE  (AMAXi.... 


^iPrsiuiit  mi 

ZERO-LIFT  CRAG  cQEEFIC 1  ENT < 
CRITICAL  MACH  NC  (FOR  CL  - 


22.0000 
2.25000 
2.00C00 
10612.0 
738.000 
17. 0000 
2.00C00 

1000.00 

1000.00 

0.002309 

1113.04 

-.174500 

7.00000 

1.40000 

•fctlSlS 

0.  QiQ750 
0.  1*0000 


INITIAL  FORWARD  VELOCITY  (  KTl  •  •  ••  120. COO 


COEFFICIENT  CF  THRUST, 

INOUCED  VELOCITY . 

DISC  AREA . . 


ti  Floss’ factor "!!! 

GROUND  EFFECT  RATIO, 


•  554153E-02 
38.8469 
1522.62 

•  650197  £-01 
•94?3o2 
1.00000 


***** ********* ******* *********************  ******* 
FORWARD  VELOCITY  IN  KNOTS  *  120.000 

******** ************************* **************** 


PARASITE  DRAG  *  806.258 

DISK  PLANE  ANGLE  CF  ATTACK 


-.759760E-01 


DYNAMIC  PARAMETERS 

All .  0.625400 

A12 . 0.806926 

T 1 .  0.467602 

T  2 .  0.31914 2 

T3 . 0.216296 

T  4 . 0.125826 

INFLOW  RATIO  ( LAHBCAi  ............  -0.040381 


stall  pgwer  calculations 

. . • . -.iC6778E-0i 

INBOARD  STALL  CORRECTION  FACTOR*.....  .0 
STALL  POWER  COEFFICIENT...... . 0 


ANGLE  OF  ATTACK  CALCULATIONS 

LONGITUDINAL  CYCLIC  ANGLE .  -0.109613 

LONGITUDINAL  CCLLECTIVE  ANGLE...  0.314954 

ALPHA  I  90  )  JOES). .  -0.048066 

ALPHA  ( 270 )  (DEC) .  12.511994 


HIGH  SPEEC  MACH  EFFECTS 

ADVANCING  ELAOE  TIP  MACH  NUMBER...  .8451 35 

CRITICAL  MACH  NUMBER.  •  .  72054  3 

DRAG  DIVERGENCE  MACH  NUMBER . 64591SE-01 

COMPRESSIBILITY  POWER  COE Ff . 521491E-Q4 

****»•*** ************* ******************************** 

POWER  REQUIRED 

***********  *******************************  ************ 


INDUCED  PC*ER  « 
PROFILE  FC.ER  « 
PARASITE  POWER  ■ 

COMP RESIBILITY  POWER 
STALL  POWER  » 


151.545 

290.464 

297.193 

•  54  *  00 1 

•  0 


TOTAL  POWER  REQUIRED  *  879.112 

************************************************* 

FORWARD  VELOCITY  IN  KNOTS  *  130.000 

****«****«  **** 4**«  ******************************* 


PARASITE  DRAG  *  94*. 234 

DISK  PLANE  ANGLE  CF  ATTACK 


*  -»8916o4E-01 


DYNAMIC  PARAMETERS 

All..... . . . . .  0.680172 

A 12  ••  . . .  0.8808  77 

. .  0.625882 

A 14..  •••••••  •••••••••  .•••••••••«•  1.20  74/4 

Tl.« . •••••• .  0.470875 

T  2.  . . . . . .  0.  325344 

T  3 . . .  0*221234 

T  4 . . .  0.136798 

INFLOW  RATIO  (LAMBDA) .  -0.  04*6  6  8 
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STALL  POWER  CALCULATIONS 

. . . . . . .  7984162-02 

INBOARD  STALL  CORRECT  ICN  FACTGR . 0 

STALL  POWER  COEFFICIENT.. . 0 

ANGLE  OF  ATTACK  CALCULATIONS 

LONGITUDINAL  CYCLIC  ANGLE . -0.122350 

LONGITUDINAL  COLLECTIVE  ANGLE...  0.326149 

ALPHA  (90  )  (OEG  )••  •  . .  -0.291443 

A  LPHA  (  27  0 1  (DEC)..... .  13.728103 


HIGH  SPEEC  MACH  EFFECTS 

ADVANCING  eLAOE  TIP  MACH  NUMBER.. •  .860309 

CRITICAL  MACH  AUMBER . .  .723294 

DRAG  DIVERGeNC  E  MACH  NUMBER . 770157E-01 

COMPRESSleUITY  POWER  COEFF.... . 630606E-U4 

it***************************************************** 

POWER  REQUIRED 

********************** ******************************** 


INDUCED  PC*ER  »  139.900 

m.m 

^MPRES^SJLITV  POWER  «  162.039 

TOTAL  POmER  REQUIRED  «  988.636 


************************************************* 
FORWARD  VELOCITY  IN  KNOTS  »  140.000 

********************* ********************* ******* 


PARASITE  DRAG  »  1097.41 

DISK  PLANE  ANGLE  OF  ATTACK  *  -.103412 


DYNAMIC  PARAMETERS 

All . 0.735630 

A 12.  •  ......  . . . .  0.  956563 

A13..  •  0.679659 

A 14 . 1.242629 

T 1 .  0.474411 

T  2 . 0.332042 

T3 . 0.224407 

T  4.  . .  0.147888 

INFLOW  RATIO  ( LAM8CA) . -0.050098 
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STALL  PCM  Eft  CALCULATIONS 


3897Q1E-02 


RT . . 

INBOARD  STALL  CORRECTION  FACTOR . .  .0 

STALL  POWER  COEFFICIENT . . . 0 


ANGLE  OF  ATTACK  CALCULATIONS 

LONG1TU01NAL  CYCLIC  ANGLE .  -0.136576 

LONGITUDINAL  COLLECTIVE  ANGLE...  0.339933 

ALPHA  (90  I  (CEG) . . . -0.520492 

ALPHA  12701  (DEG) . 15.129103 


HIGH  SPEEC  MACH  EFFECTS 

ADVANCING  ElADE  TIP  MACH  NUMBER...  .875484 

CRITICAL  MACh  NUMBER.  ...  . 725882 

ORAG  DIVERGENCE  MACH  NUMBER . 896012E-G1 

COMPRESS  I  El LIT Y  POKER  COEFF . 745875E-04 

*************** ************************** ************* 

PCWER  REQUIRED 

******************* ***  ******************************** 

INOUCED  PCkER  «  129.913 

PROFILE  PCwER  *  322.449 

PARASITE  FCWcR  *  471.789 

COMPRESIBILITY  POKER  »  191.657 

STALL  POWER  »  .0 

TOTAL  POWER  REQUIRED  »  1115.81 

**************  *********************************** 
FORWARC  VELOCITY  IN  KNOTS  »  150.  OCO 

*  *************  *******  *************************  *** 

PARASITE  DRAG  «  1259.78 

DISK  PLANE  ANGLE  OF  ATTACK  «  -.118713 


DYNAMIC  PARAMETERS 

All . ••••••••• . .  0.791849 

A  12 . . . .  1.  03*1  7  C 

A 13 .  0.754801 

A 24*.  •  ••  •  .  1.  281052 

. . ••••  0.478208 

T2 . •••••••  0.339237 

T3 .  0.227815 

T4.  ••  •  ••  . .  0.159103 

INFLOW  RATIO  (LAMBDA)..... . -0.056720 
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I 


STALL  power  calculations 


AT 


•  2307  77E-02 


GAMMA 

if 


VALUES  FOP  DETERMINING  STALL  COEFFICIENT 


-C.264593S0 
-0.14755$; 
0.32< 


>449< 


INBOARD  STALL  CORRECTION  FACTOR . 662213C94 

STALL  POwEP  COEFFICIENT . 312025956-3* 


ANGLE  OF  ATTACK  CALCULATIONS 


CYCLIC  ANGLE .  -0.152545 

_  COLLECTIVE  ANGLE...  0.356377 

ALPHA  (  90)  (CEGJ . . .  -0.738637 

ALPHA  (270J  (DEG).  •  ...  . . .  16.740753 


LONGITUDINAL 

LONGITUDINAL 


HIGH  SPEED  MACH  EFFECTS 

ADVANCING  eLADE  TIP  MACH  NUMBER...  .890o58 

CRITICAL  MACH  NUMBER . 728348 

DRAG  DIVERGENCE  MACH  NUMBER . 102309 

COMPRESSIBILITY  POWER  COEFF. .  . . 867333E-04 


****************************************************** 

POWER  RE  GUI  RED 

****************************************************** 

INDUCED  PChER  »  121.281 

PROFILE  PChER  «  336.941 

PARASITE  POWER  »  580.279 

CQMPRESI  BIHT  Y  POWER  »  223.  009 

STALL  POWER  >  80.1  774 

TOTAL  POWER  REQUIRED  *  1341.69 


**************  *********************************** 
FORWARD  VELOCITY  IN  KNOTS  «  160.000 

**************  *******  **************************** 


PARASITE  ORAG  *  1453.35 

DISK  PLANE  ANGLE  OF  ATTACK  >  -.135063 
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DYNAMIC  PARAMETERS 


All . .  0.  8489 C 7 

A  12 . ..  1.113397 

A  13 .  3.791449 

A 14 . 1.312859 

Tl.... . 3.432267 

T  2 . 3.346927 

T3 .  3. 2214: 7 

T  . . 3.170453 

INFLOW  RATIO  (  LAMflOAJ . -0.064591 


STALL  POWER  CALCULATIONS 


•  1 15352E-01 


VALUES  FOP  DETERMINING  STALL  COEFFICIENT 


GAMMA 

CS 

BS 


-C.  30194  568 
-C.  17515445 
0.36584556 

0.7398598 

1.3566751 


INBOARD  STALL  CORRECTION 
STALL  POWER  COEFFICIENT.. 


FACTOR .  1.30000  000 

•••..•.••a..  .635492* Cc-34 


ANGLE 


ATTACK  CALCULATIONS 


LONGITUDINAL  CYCLIC  ANGLE .  -3.17C535 

LONGITUDINAL  COLLECTIVE  ANGLE...  3.375743 

ALPHA  (90)  ( JEG  ) .  -0.949361 

ALPHA  ( 270  )  (CEO . 18.591431 


HIGH  SPEEC  MACH  EFFECTS 

ADVANCING  6 LAO E  TIP  MACH  NUMBER...  .905332 

CRITICAL  MACH  NUMBER.  . . . 73C729 

DRAG  DIVERGENCE  MACH  NUMBER . 115133 

COMPRESSIBILITY  POWER  COEFF . 997223E-04 

****************************************************** 

POWER  REQUIRED 

***************  *************************************** 


INDUCED  POWER  * 
PROFILE  PCnER  * 
PARASITE  POWER  * 

COMP  RES  I  B  HIT  Y  POWER 
STALL  POWER  » 


113.722 
352.433 
7 04.  243 
25o . 243 
155. 586 


TOTAL  POwER  REQUIRED  «  1582.23 
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***** ***************************************** *** 

FCRaARO  VELOCITY  IN  KNOTS  *  163.359 

**************  ******* ******************** ******** 


PARASITE  DRAG  ■  14*4.16 

DISK  PLANE  ANGLE  OF  ATTACK  *  -.140799 


DYNAMIC  PARAMETERS 

All . . .  0.  863277 

. . 1.141191 

A 13  .  .  •  •  . . 0.810842 

Alt .  1.337756 

T 1 . . .  0. *83685 

T2*. . 0«349o22 

. .  0.232734 

T4.  ••  •  . . 0.174297 

INFLOW  RAT1C  <  LAMSCA) . -0.067526 


STALL  POWER  CALCULATIONS 
RT . •••• .  15598  7t-0 1 


GAMMA 
CS 
8  S 

xs 

xo 


VALUES  POP  DETERMINING  STALL  COEFFICIENT 

-C.31570089 
-C. 18555188 

s  n.  3  £06 -s  A?  ■ 


0.  38093823 

0.  733  6488 
1. *493  771 


INBOARD  STALL  CORRECTION  FACTOR .  1.00000000 

STALL  POWER  COEFFICIENT . 611C876CE-04 


ANGLE  OF  ATTACK  CALCULATIONS 

LONGITUDINAL  CYCLIC  ANGLE . -0.177082 

LONGITUDINAL  COLLECTIVE  ANGLE...  0.3829*7 

ALPHA  ( 90  )  ICEG  I. ........  .. .  -1.019017 

A  LPHA  1270)  I  DEC-). . 19.272018 


HIGH  SPEEC  MACH  EFFECTS 

ADVANCING  8LA0E  TIP  MACH  NUMBER...  .910929 

CRITICAL  MACh  MJMBER . . . 731516 

DRAG  DIVERGENCE  MACH  NUMBER. ••••• •  .119*12 

COMPRESSIBILITY  POWER  COEFF . .  .1042*1E-0: 
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****  *******  4***  *******************************  *******  V 


PCaER  REQUIRED 

****************************************************** 


PARASITE 

hew 


PCnER  * 

PCWcR  * 

_  POWER  * 

_ 28ILZT Y  POWER 
POWER  « 


HfcHf 

749.536 
267. S55 
157.023 


TOTAL  POWER  REQUIRED  »  1643.66 


INITIAL  STALL  CMS E T  VELOCITY  APPROXIMATION  UTSJ  «  76. 


DIFFERENCE  BETWEEN  A270  AND  AMAX  *  .259871 E-02 

ACCEP4TA8LE  ERROR  *  .10  0  0  00E-02 


**************  *******  ************************** 


ANuLE  AT  STALL  ONSET-  13.9811 

STALL  CNSET  VELOCITY  -  131.913 

VELOCITY  MAX  Ef^DURANC  E  »  76.0347 

MAXIMUM  FCR-ARC  VELOCITY  -  io3.359 

*****************************************  **<•*** 
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G.  COMPARISON  OF  PROGRAM  OUTPUT  VS  TEST  D AT  A 

This  section  compares  the  output  cf  the  HP41-CV  ar.c  IBS 
3033  computer  programs  to  actual  flight  test  data  gathered 
at  the  Naval  Weapons  System  Canter. 


POWER  REQUIRED  (HP) 

1*0.0  ***.*  7*#.#  i***.*  in*.*  1000.0  Ira*.* 


POWER  REQUIRED  VS  VELOCITY 
FOR  THE  AH1-J  USING  THE  HP  41-CV 


Figure  c.1.  Poeer  Carres  Senerat9d  by  HP*1-CV  Progr 


POWER  REQUIRED  (HP) 

t.i  t«o. •  •it.i  i  loii.i  itu.o  mt.a 
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POWER  REQUIRED  VS  VELOCITY 


Pigure  C.  3, 


Comparison  Between  Computer  Data  and  Actual 
Test  Plight  Data. 
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